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O ■ Abstract: Motivated by a heuristic model of the Yang-Mills vacuum that ac- 

curately describes the string-tension in three dimensions we develop a systematic 

J-.- method for solving the functional Schrodinger equation in a derivative expansion. 

^ ■ This is applied to the Landau- Ginzburg theory that describes surface critical scaling 

in the Ising model. A Renormalisation Group analysis of the solution yields the value 



7] = 1.003 for the anomalous dimension of the correlation function of surface spins 
which compares well with the exact result of unity implied by Onsager's solution. 
We give the expansion of the corresponding /5- function to 17-th order (which receives 
contributions from up to 17-loops in conventional perturbation theory). 
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1. Introduction 

In recent years there have been considerable advances in our understanding of the 
quantum field theory of systems with higher conserved charges. Deep mathematical 
structures have emerged in the study of integrable models in two dimensions and in 
supersymmetric theories in higher dimensions. However the field theories that are 
directly applicable to the physical world do not manifest the high degree of symmetry 
needed for mathematical tractability. This limits us to the use of perturbation the- 
ory, which cannot describe strongly interacting theories, and the largely numerical 
approach of lattice field theory. This suggests that it is worthwhile to develop new 
methods for the solution of generic non-integrable quantum field theories. In this 
paper we will describe an approach to quantum field theory in which the functional 
Schrodinger equation is constructed in an essentially standard way, [|^-[]TB[, but is 



then solved using a derivative expansion for the wave-functional from which the con- 
tinuum limit is obtained using a version of the Borel transform. Like perturbation 
theory this yields an approximation in the form of a series, but unlike perturba- 
tion theory the method applies at arbitrary values of the coupling. The approach is 
inspired by a heuristic model of the Yang-Mills vacuum but we will develop it for 
the toy-model of scalar ip^ theory in 1 + 1 dimensions. This is partly because the 
simplicity of the model allows the series expansion to be calculated to reasonably 
high order using modest computing resources but also because the strong-coupling 
behaviour of the vacuum wave-functional describes the surface critical scaling of the 
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Ising model and so we are able to test the method against exact results which derive 
from the Onsager solution and conformal field theory. 

In a D-dimensional gauge theory there is a constant force between static quarks 
leading to confinement when the Wilson loop, iyL[C], associated with a closed curve, 
C, depends on the minimal area spanning C as exp(— cr Area [C]), |]T^. a is the string 
tension. WL[C] is the vacuum expectation value of the trace of the path-ordered 
exponential of the gauge-field, WL[C] = (0|trP exp{— J^ A ■ dx) |0). This has a 
functional integral representation J 'DA exp (— 515/(7^) trP exp(— J^A ■ dx) where 
Sd is the Yang- Mills action in D Euclidean space-time dimensions, and g a coupling. 
Areal behaviour is observed in Monte Carlo lattice studies for D = 2, 3 and 4, but 
has not been obtained by analytic calculation except in two dimensions where Yang- 



Mills theory becomes free with the gauge choice Ai = 0. In 1979 Greensite [18 



and Halpern [|T9| proposed a simple model of the confining Yang-Mills vacuum. For 



simplicity consider D = 3. If we choose the temporal gauge Aq = and work in 
the Schrodinger representation so that Ai and A2 are diagonal on the quantisation 
surface t = then the representative of the vacuum {A\0) = "^[A] is a gauge- invariant 
functional of 74i(a;^,x^) and A2{x^,x'^). The idea is to model this three-dimensional 
vacuum in terms of the two-dimensional action as \E' ~ exp{—kS2/g'^), with k a 
dimensionless constant. This guarantees gauge-invariance and if we choose a planar 
C so that it can be rotated into the quantisation surface then the usual expression 
of vacuum expectation values within the Schrodinger representation yields 

WL[C]= f PAiDA2^*[v4]trPexp(- f A-dx)'^[A] 

= VAiVA2exp{~2kS2/g^)tTP exp{- A-dx) (1.1) 

which reduces to the Wilson-loop for a two-dimensional gauge theory for which we 
can compute the string tension in terms of k. This form of ^[^4] was compared 
directly to Monte Carlo simulations of the vacuum functional in lattice gauge theory 
in [^. It might now be timely to repeat this comparison given the improvements in 
lattice gauge theory technology. This argument was generalised to quantum gravity 

in p-m. 

"^[A] should satisfy the Schrodinger equation, and this has been used to estimate 
k, [p. The Hamiltonian is — / d'^xtr {A^ + B^) where the 'magnetic field', B, is the 



1, 2-component of the field strength from which the two-dimensional action is con- 
structed as 5*2 = — J d'^xti {B"^). In the Schrodinger representation A = —ig'^6/6A, 
so that the Schrodinger equation is 

H^ = f d^xti L^—^-^bA ^ = E^. (1.2) 

Differentiating S2 functionally with respect to A gives the covariant derivative of B, 
DiB. Differentiating again gives a second order partial differential operator, Q, so 
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H^ = [ ^Tin + ^k'^ / d^xtT{{DBy) + ^S2 ) ^ = E^. (1.3) 



that for the model vacuum the Schrodinger equation becomes 

-^Trfi + ^A;^ f d\ ii {{D Bf) + - 

y 9 J 9 

When the functional trace of Vt is regulated using a heat-kernel based regulator to 
preserve gauge invariance the result is a divergent constant plus a multiple of 5*2, so 
that fc, and hence a, are determined if we ignore (DB)"^ as a higher order correction. 
Specifically, if zeta-function regularisation [^ is used we obtain 



1927r ^ ^ ^ 



for the gauge-group SU{N), (and where we have adapted the result of [^ to the 
convention tr [TaTi,) = —5ab/'^ for the generators of the fundamental representation.) 
This formula shows remarkable agreement with recent lattice studies. Figure (|l]) plots 



^/o^/{g'^N) against A^ as measured in Monte Carlo simulations [Q. ^/a^/{g'^N) is 
shown on the same figure, and the values agree to within a half to one percent which 
is significant given that the formula ( |1.4| ) has no adjustable parameters. (There is a 
similar formula, axKN = 9'^{N'^ — l)/(87r), |^, based on a different parametrisation 



of the canonical variables which is also displayed in the Figure). 

This model begs many questions: Can the vacuum functional really be the ex- 
ponential of the integral of a local quantity, i.e. depending on a single point in space, 
since this is quite unlike the non-local functionals seen in perturbation theory. For 
example the free Maxwell field in four- dimensions has 

^ = exp(- / d?x(fyB{x) B{y)/{e^\x - yH) . 

Is the approach regularisation dependent, and is it consistent with the usual ultra- 
violet behaviour and renormalisation of gauge theories? Is the model Lorentz in- 
variant? Can corrections to the model be computed and shown to converge to the 
lowest energy solution to the Schrodinger equation? Applying the argument to four- 
dimensional gauge theory would give a vacuum functional that is the exponential 
of the three-dimensional action with a cut-off dependent coefficient that appears to 
lead to a divergent string tension. Despite these issues the agreement between lattice 
simulations and the formula for the string tension suggests that it is worthwhile to 
investigate the approach further and make it more systematic. 

In calculating k the contribution of {DB)'^/g^ was ignored in ( |1.3|) . This may 
be corrected for by including a similar term in the vacuum functional. This leads to 
an adjustment to the coefficient of the original term proportional to tr(i?^) because 
J dPxd"^ /5A^ applied to / dPxti^^DBY) gives a piece proportional to e^^ / (Pxti B"^, 
where e is a short-distance cut-off as well as j (fxtii^i^DBY) so that now, if we set 
\E' = exp W then the approximation for W is 
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Figure 1: ^/{g^^/N"^ — 1) against N from lattice QCD, ||2^. The horizontal hnes 
correspond to the predictions from o"^ (lower) and ctkkn- 



which satisfies the Schrodinger equation for appropriate values of ^2 and k^ if now 
we ignore terms of mass dimension eight. These higher dimension expressions can 
then be taken into account by adding further terms so that the process generates an 
expansion for W consisting of local expressions of increasing dimension with coeffi- 
cients that are themselves series in decreasing powers of the cut-off e. Superficially 
such a derivative expansion would appear to be applicable only to slowly varying 
fields, making it difficult to address the ultra-violet behaviour that has to be un- 
derstood in order to renormalise the theory. Furthermore the cut-off dependence of 
the coefficients is unexpected since a theorem due to Symanzik [^ implies that W 
should be finite as e — ;> (at least to all orders in perturbation theory). 

To address the problem of whether a systematic expansion of the vacuum func- 
tional can be obtained as a derivative expansion with the correct ultra-violet be- 
haviour we turn to a simpler model than gauge theory, namely scalar (p^ theory in 
two space-time dimensions. The simplicity of the model will allow us to pursue the 
expansion to high order and, as stated earlier, the strong-coupling ultra-violet be- 
haviour is related to the surface critical scaling of the Ising model, providing exact 
results to test our approach against. 
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The Hamiltonian for scalar ip in two space-time dimensions is 

H = jdx{^- {^l + ^'^ + ml ^g) + I ^^) , (1.5) 

where (^o, ""^o and qq are the bare field, mass and coupling. Quantisation consists 
simply of replacing 0o by —i6/6ipQ. We can now look for eigenfunctionals of H. We 
will assume that, as in quantum mechanics, the vacuum functional has no nodes, so 
that it can be written as \E' = exp W, giving a Schrodinger equation 



^^\-^^- \^] +'^0 +^o¥>o + T7;'^o] =Ev (1.6) 



2 

S^l \6ipo) ' ''' ' "^'"'' ' 12 

where -E^, is the energy of the vacuum. The first term needs to be regularised, and 
we adopt a simple momentum cut-off, replacing S^/Sipl by 

/, = fdxdy f ^±e^K^-y) f (1.7) 

with s a small distance. 

The problem of renormalisation is to tune the bare quantities in terms of s to 
obtain finite physics. Renormalisation conditions have also to be chosen to replace 
the (potentially divergent) bare quantities by finite parameters, in this case a mass 
scale /i and dimensionless coupling g. We are aided in this by Symanzik's theorem 
pO|] which implies that provided we renormalise ipQ multiplicatively Lp^ = vZ (f (and 
tune mo and f^o) then W is finite apart from a local term proportional to J dxip"^. 

Conventional semi-classical perturbation theory can be developed |jl| by re- 
instating h in the quantisation 00 = —ih6/6ipo and expanding W = Wi/h + W2 + 
hWs + .. so that to leading order ( p.. 6] ) becomes the Hamilton- Jacobi equation with 
a solution most conveniently expressed in terms of the Fourier transform if as 

Wi = - —<f{-p)uj{p)<f{p)-— / -—..--27c6{pi + ..+p4) +.. 

y 47r 4! y 27r 27r uj{pi) + .. + uj{pi) 

(1-8) 

where u = \Jp^ -|- mg , and bare and renormalised quantities coincide at this (tree) 

level. Clearly for slowly varying ip we can take \p\ << itlq and expand uj = tuq + 

p^/(2?Tio) + •• to obtain a local expansion. The one-loop, i.e. order h, correction 

can be readily obtained, ||l|, |^, but this is not our purpose here. Instead we will 



develop the approach of [^ and consider solving for the full W starting from an 
expansion in a basis of local functionals of the renormalised field with coefficients 
{b}: 



W 



I dx (62<^' + &4<^^ + &6<^^ + •• + &2,2<^" + &4,2<^'V + - + &2,4<^"' + ••) (1-9) 



-5- 



Now the regulated second functional derivative acts on this local expansion to pro- 
duce another local expansion, for example 



dx 



^' ^ dy^\'^= [dx fAy,'2^ 2 ^2| ^^^^Q^ 



S(^' 



\ ITS Svrs 



The other terms in the Schrodinger equation are also local functionals of (p so we can 
reduce this functional differential equation to an infinite set of algebraic equations for 
the infinite number of constants {b} by equating coefficients of the basis functionals 
to zero. To solve these requires some form of truncation, for example by simply 
ignoring coefficients of local terms containing more than a prescribed number of fields 
or derivatives. In [^ a different truncation was proposed that would use an estimate 



of such 'high order' coefficients. Once the problem has been made finite the equations 
can be solved, for example numerically, although the stability and convergence of the 
solutions may be sensitive to the truncation. We will use a different approach to that 



of |3^ . Rather than truncate the expansion ( |1.9| ) we look for solutions to the full set 
of (untruncated) equations in the form of expansions of the coefficients {b} as power 
series in 1/s. We will see that the enlarged set of equations obtained by equating 
to zero the various powers of s form closed systems that can be solved rapidly and 
simply in terms of the renormalised quantities /x and g. The quantities Z, g^ and tuq 
will similarly be obtained as power series in 1/s. We will compute these series up to 
some fixed order, N, determined by computing resources. This is the only truncation 
we will use. 

The continuum limit is obtained by sending s to zero, (and A^ to infinity) so 
these series solutions would appear to have little physical relevance, as they are only 
likely to converge for large s. However the small s behaviour can be extracted using 
a version of the Borel transform, which we describe in Section 3. 

We will test our solution of the Schrodinger equation against exact results ob- 
tained for the surface critical scaling of the Ising model from the Onsager solution 



and conformal field theory |^T|. At a critical value of the renormalised coupling, g, 



the /i — > limit of cp^ theory describes the Ising model at criticality. More specifically 
the wave-functional describes the Ising model with a boundary (corresponding to the 
quantisation surface) at the 'Ordinary Transition'. That part of the wave-functional 
that is quadratic in p is related to the correlator of two spins lying in the boundary 
which is known to depend on their separation |x| as l/|a;|. This quadratic piece 
can be written as J dx dy ip{x)(f{y) /\x — yl"^'^ where t] is the anomalous dimension, 
so the correspondence with the Ising model requires that r] = 1. Now a standard 
Renormalisation Group argument presented in Section 4 gives rj as the value of 
(i log Z/(i log /i = 7(5') computed at the zero of the beta-function dg/dlogfi = (3{g). 
In Figure ^ we display these two functions computed from our solution when we trun- 
cate the series expansions to order 1/s^^ terms. The graph indicates good agreement 



-6- 



1- 
























-1- 


■ 






















-•- 




- 




















+ 






- 


















+ 








. 


















+ 








0.5- 














-•- 


+ 










■ 












+ 














- 










-1- 














-♦■ 


. 








+ 


















■ 


-1- 


+ 


0.5 






1 






1 .5 






2 


o~ 


























■ 






















■•- 




■ 


-1- 
























0.5 - 


























■ 




-•- 
















-f- 






. 






-1- 




















- 


















-1- 








— 1- 








+ 


















- 










+ 








-•- 








- 












+ 


-1- 


+ 











Figure 2: ^{g) (upper curve) and I3{g) from the Schrodinger equation. 

with the exact result, and as we will see the agreement is considerably improved 
upon extrapolation in the number of terms, N. 

2. Solving the Schrodinger Equation for cp^. 

The Schrodinger equation for the vacuum, HexpW = E^explV, can be written in 
terms of the renormalised field as 



dx 



6^W 



5ip^ 



+ 



6W 

6ip 



- ZV" - m^zV' - — ^V' + 2Z£. ) = (2.1) 



where S^ is the energy density of the vacuum. The local expansion ( |1.9| ) takes the 
form 



W 



I dx( 62 <^^ + &4 V^ + &6 <^^ + •• + h,2 <^'^ + &2,4 '^"^ + •• 

+&4,2 V^'^ ^"^ + &4,4,1 V^"^ V'^ + &4,4,2 ^"^ + ■■ 

+be,2 <^'V + •• + h <^' + &8,2 <^" <^' + ..) (2.2) 



where the dots denote an infinite number of terms containing increasing powers of 
the field and its derivatives. The two (or three) labels on the coefficients denote the 
number of fields, the number of derivatives (and a further label to distinguish between 
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different terms for which the first two labels are the same). We will not truncate the 
expansion at this stage. Substituting (|L^) into the Schrodinger equation gives 



dx 



262 262 2 
ITS Stts^ 



+ 



,2 I 2&4.2 , 464 4.1 OC4 4 2 , \ , 4 

' H 4- + ^ o + ■■] + ^ 



8bA 



1264 



TTS 



+ 



26. 



4,2 



Svrs" 



ITS 

5668 



TTS-' 



Svrs'^ 



+ .. 



3066 ^ 266,2 ^ 



vrs 



+ 



2^8,2 , A , 2/2 /^ 1265,2 






¥^^462 + ip'^8b2b2,2 + ¥^^166264 + y^'^ (246265 + 166^) + 
(/?V'^ (166264,2 + 4862,264) + .. 



r72 /2 2 ry2 2 90 ry3 4 , r, ry c 

-Z (f — ttIqZ (f Z ip + 2Zty 







(2.3) 



We can also look for an eigenstate of the Hamiltonian describing a single particle 
at rest in the form \E' = [/ exp W The Schrodinger equation for this state, H'^ = 
{fi + Ey)"^ , can be written as 



dx 



6^U 



Sifi^ 



+ 2^f + 2,Zt/Uo 

dip dip J 



(2.4) 



Take U as another local expansion: 

V = dx(ip + Csp^ + C^ip^ + .. + C3^2^p''^ + CsA^p"^ + .. + C5,2(/?V^ + ••) > (2-5) 

where again the dots denote an infinite number of local terms all of which we keep in 
the expansion for the time being. Substituting into the Schrodinger equation gives 



dx 



Qc^ 



2c, 



2c, 



3 ^L^3,2 ^'-3,4 I \ I 3 



I 12 I 6C5,2 



2OC5 2C5,2 



TTS 



Svrs'^ 



TTS 



+2 



(^262 + p (64 + 6C362) + v^v?' (12c362,2 + 6c3,262 + 264,2 



+2/iZ (y? + C3(/?^ + C5P^ + .. + C3,2¥^¥^'^ + C3,4¥^(/'"^ + ..))= (2.6) 



These equations have to be solved in conjunction with renormalisation conditions 
that identify parameters that remain finite as the cutoff is removed. We will take /i 
to be our finite mass-scale. Symanzik's theorem [^] implies that the coefficients 6 in 
the local expansion of W are finite, with the possible exception of 62, so we choose 
renormalisation conditions 



h 9f^ . 1 



[2.7) 



with dimensionless g playing the role of renorinalised coupling. Our aim is to solve 
( p.3|) and (|2.6|) for the coefficients h and c, and also for qq^ itlq and Z in terms of the 



finite parameters /x and g. 

Equating the coefficients of the basis functionals to zero in (2^) gives an infinite 
set of algebraic equations. This may be divided into two classes. In the first class 
are the equations obtained from the coefficients of 1, y?^, v^'^, and ip^. 

262 2622 ^ 12^4 2642 ,999 . ^ 

TTS 6TTS-^ ITS SlTS-^ 

\ ^ + ^ ^ + .. + 86202,2-^ =0, hr ^ + .. + 166204-777^ =0 

(2.9) 
These equations involve the bare quantities go, mo, Z and S^ as well as the coefficients 
b. The second class consists of all the other equations, and these involve only the 
coefficients b. 

^66 26 126 

— - + T^^ + ■■ + 246266 + 166^ = 0, — ^ + .. + 166264 2 + 4862 264 = 0. (2.10) 

ITS Stts-^ its 

(2.11) 

We begin by solving this second class for the 6 coefficients other than 62, 64 and 
62^2 in terms of g, fi and 62 using ( p.7| ). Expanding each coefficient as a power series, 
bn = ^r=o 00 K/^^y ^^^ then equating the coefficients of powers of 1/s to zero in 
each equation gives further equations that can be easily solved in MAPLE to high 
order. An important feature is that they form a closed system which is easily solved 
order by order in 1/s. Consider the equation corresponding to the coefficient of a 
given local term involving F powers of the field and a total of D derivatives. The 
equation obtained by equating to zero the coefficient of s~^ takes the form: 



UP 



D+p-1 

E^P,D',r' jr,~l+D~D' , > >P 
^F,D,r "F+2,D',r' + ''S 0^ 

D'=D 

D"=D,F'=F-4,p'=p 
, "SP ^p,D",F',p',r',r" jv' jv-p' _ r. /^ , „n^ 

-r 2_^ ^F,D,r "F',D",r' "F+2-F',D-D",r" ~ '^ l^-J^^J 

D"=0,F'=4,p'=0,r',r" 

where A and C are known numerical coefficients. This can be solved for ^^^ in 
terms of the coefficients {6^,_^,^,} where F' + 2D' + 2p' < F + 2D + 2p = Q and 
either F' = F,p' <p or F' < F,p' < p, D' < D. 

D+p-l 



LP L/^ V^ Ap,D'y jv-l+D-D' 

^F,D,r — 7 y Z_^ ^F,D,r '^F+2,D',r' 



^2 D'=D 
D"=D,F'=F~A,p'=p 

Sr^ ^T) n" F' -n' r' r" 



. V^ ^p,D",F',p',r',r" jv' jv-p' \ ^o i oN 

' /_, ^F,D,r "F',D",r' "F+2-F',D-D",r" y^-i^o) 



D"=0,F'=4,p'=0,r',r" 
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To obtain fe^^Jr explicitly we can build up the solutions starting with p = 0. Given 
that the renormalisation conditions fix 6° and 622 ^^ begin by calculating 6g, 634 
and 642 and then use these results to calculate 6g, 644, 6g2 ^^^ fe^e- Fro^i these we 
compute more coefficients until all the {fe^^, ^,} are known with F' + 2D' < Q. Then 
we move on to {bp, ^i} by first using the results obtained to compute 6g, 634 ^^^ 
b\ 2- From these we obtain 6g, 64 4, 6g 2 and b^ g and continue the process until all the 
{b]p, £,,} are known with F' + 2D' + 2 < Q, and from these we can compute {6|,, ^,} 
and so on. Although this appears complicated in practice it is straightforward. For 
example ( |2.10|) gives to lowest order in 1/s 



24626° + ^^ = 0^6° = -M!, 16626^2 + ^ = ^ 6°2 = -^. (2.14) 
^4 ^ 9662 ^ ^'2 2 ^'^ 3262 ^ ' 

Once the {&^£)} are known substituting the solutions into the first class of equa- 
tions yields series for the bare quantities. Thus the first of ( p.9| ) gives 

n 262 3g 

Z' = --T7^ + - 2.15 

/i lo627rs 



Having obtained the 6-coefficients and Z, tuq and go we solve ( p.6|) to obtain the 
c-coefficients and the energy eigenvalue in terms of /i, g and 62 again as power series 
in 1/s. Finally, given that the energy eigenvalue equals fi we can invert the series to 
obtain 62 in terms of g and n, and substitute this into all the other series expansions 
to express all quantities as power series in 1/s depending on g and /x. In the Appendix 
we give the expansions ^ for g^, and Z up to order 1/s^^ which requires taking terms 



with F + 2D up to 40 in the local expansions ( p.2|) and ( p.5|) , of which there are 
about 625 and 540 respectively. 

These expansions can only be expected to converge for large values of s, if at all, 
but the continuum limit requires taking s to zero. In the next section we will show 
how a version of the Borel transform can be used to extract the small s behaviour 
from these series, and then in Section 4 we will apply this to the bare quantities to 
extract an anomalous scaling dimension which we can test against the exact result 
for the Ising model. 



^ We note in passing that such series could be obtained from standard Feynman diagram expanded 
in powers of 1/s, for large s. Given that these expressions contain powers of the coupling to order 
(7^^ and g^^ this would require working to 17 loops in conventional perturbation theory. 
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3. Small s from large s via Borel 

The functional integral representation of wave functional described in section 4 can 
be expanded in terms of Feynman diagrams. For massive field theories Feynman 
diagrams are analytic in the cut-off s for large enough s, although for small s they 
have singularities in the left-half-plane. After renormalisation they have finite limits 
both as s ^ and as s ^ oo (where they tend to their tree- level values). Assuming 
that these properties hold beyond perturbation theory we will exploit this analyticity 
to extract the small-s behaviour from our series expansions. 

As an example, consider a function f[s) that has the assumed analyticity proper- 
ties of wave-functionals, namely it is analytic for large s, and in the right-half-plane, 
with finite limits as s approaches the origin and infinity. For large s it will have a 
Laurent expansion containing only non-positive powers of s, but suppose that we only 
know the first few terms, An{s) = Yln=o N^n/s"", from which we wish to estimate 
the small s behaviour. We can obtain this from the integral 

where C is a large circular contour centred on the origin running from just below the 
negative real axis to just above. This integral is readily computable in terms of the 
coefficients a„. Also, since Aj^{s) ~ /(s) on C: 

If for Ix the contour C is now collapsed to surround the singularities of the integrand, 
then the pole at sq contributes /(so)- Poles and cuts lying to the left of sq will be 
suppressed when Re(so) > 0, by factors exponential in A. For example, a pole at 
s = Sp contributes a piece proportional to exp(A(sp — Sq))- Similarly the size of the 
contribution from a contour, Cf, of finite length extending as far to the right as Sc 
is bounded by 

This exponential damping means that for large A 

I^ ^ fiso) . (3.4) 

We cannot take A arbitrarily large, but only as large as allows I^ to be a good 
approximation to Ix, which is the case when I^ ^^ differs by no more than say a 
few per cent from I^ . If this condition fixes A to be A* then we arrive at the 
approximation /(sq) ~ ^x* 
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yN 



We will illustrate this with the truncated series A^ = X]n=o(~V'5)"/("'+ !)• ^s 
N —^ oo this series converges for |s| > 1 to / = slog(s/(s + 1)), which is analytic 
except for a cut from s = — 1 to the origin, so 



■50 log 



■50 



so + 1 



y 



dy- 
y + Sq 



,-Ky+so) 



(3.5) 



from which it is clear that limA^oo I\ = /(■5o)- If we use the series truncated at say 
nineteen terms then ^19(5) itself approximates f{s) to better than one percent for 
s > 1.07. At s = 1 the error grows to three and a half per cent and for |s| < 1 the error 
increases rapidly as s gets smaller. Figure |^ plots I^^ and /|^ for sq = 0.9. The two 




Figure 3: /|^ (lower curve) and /j[ 



curves agree up to about A = 9, and for this value provide a very good approximation 
to /(so) = -0.6725, (If = -0.6724 and Jf = -0.6727). For larger values the two 
curves diverge, and approach ^19(0. 9) = —0.4823 and ^13(0. 9) = —0.8524 which 
are poor approximations to /(0.9). For much smaller values of Sq we still obtain 
reasonable results, for example in the case of so = 0.05 the two curves diverge at A ~ 7 
yielding an approximation of /j[^ = —0.2123 to /(0.05) = —0.1522, which is fifty per 
cent out, but it is amusing to compare this to the value of y4i9(0.05) = 0.25 x 10^"^. 
We will be principally interested in the value of / at the origin. We can take 
So — > in the general case, but the new feature that arises is that any cut in the 
left half-plane that ends at the origin will no longer be exponentially suppressed. 
However, our assumption that /(s) has a finite limit as s ^ implies that if / ~ s'' 



-12- 



then K > and this behaviour is power law suppressed. Also the expression for I^ 
simplifies to ^„=o n in^^ /n\. In our example I\ = e~^/\ — 1/A and the second term 
demonstrates this power law suppression. Again II^ and II^ agree for A up to about 
7 yielding an approximation of —0.1416 to /(O) = 0. 




Figure 4: I^ (upper curve) and I^ for s ^ 



This approximation can be greatly improved by changing the complex variable 
s, for example by considering /(s'*) with a > 1 instead of f{s), so that now 



If = J2 a„A"7r(na + l) 



(3.6) 



n=O..N 



Increasing a reduces the size of the last term in I^ in comparison to the penultimate 
term, because of the F function in ( ^.61 ), enabling us to take a larger value of A, 
however increasing a enlarges the region of non-analyticity. Initially this will result 
in singularities occurring at values of s with large negative parts, causing oscillations 
in the plot of I^ against A. This can be seen in our example, in Figure ^. 

Maxima appear for a = 1.4, 1.5 and 1.6 due to these oscillations. Increasing a 
still further can cause the singularities to migrate to the right half-plane, spoiling 
our construction altogether. The position of these dominant singularities can be 
estimated by studying the shape of the curve, thus the exp(— A) decay in our example 
for a = 1 can be seen by fitting numerically the curve of I^ to an exponential. 
To improve our approximation we need to increase a as much as we can without 
encountering these pathologies. We can do this by taking a as large as possible whilst 
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Figure 5: /^ and /|° for a = 1 (lowest curve), 1.1, .., 1.6 



maintaining I^ as a monotonic function of A in the region where it approximates /(O), 
A < A*, (we may have to exclude a region close to the origin if Ix has a turning point 
not due to a singularity occurring at complex s, this complication can be resolved 
by studying the iV-dependence) . So in our example we would take a = 1.348004 and 
obtain the approximation 0.013 when we require that I^^ and I^^ differ by 0.01 at 
A* = 9.77. 

The precision with which we require ll^ and I^f" to agree introduces arbitrariness 
into A*, however because we tune a so that /j[^ is quite flat for A ~ A* the estimate 
we obtain is fairly insensitive to this choice of precision, so for example if we require 
that ll^ and I^^ differ by only 0.001 (instead of 0.01) at A* the estimate changes to 
0.0115 even though A* has changed to 8.93. Note that /j[^ is considerably less flat 
than ll^ in the vicinity of A*. It has a single maximum in this region, whereas the 
flat portion of the latter turns into two turning points that are close together if a is 
slightly increased. To minimise arbitrariness we should in general base our estimates 
on the flatter curve, which means using I^ with A^ odd in this example. 

The integral Ix is the Borel transform of the function / centred at sq. Our 
approach based on studying the large-A behaviour to reconstruct the original function 
has the advantage over the usual method of avoiding the need to analytically continue 
the transform, (for example using Pade approximants), prior to inversion using the 
Laplace transform. 

In our application to the Schrodinger equation s is a sharp momentum cut-off. 
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Figure 6: I^ (upper curve) and ij^^ for a = 1.348004 in the vicinity of A^ 



It has been observed that in the different, but related context, of Polchinski's flow 
equation that a smooth cut-off increases the rate of convergence of a series solution. 
The reverse occurs here. Suppose that the series Aj^ resulted from an integral with 
a sharp cut-off as might be the case for example for a one-loop Feynman diagram. 
Thus 

N i/s N 

An = Y^ (^nls^ = rffc ^ a„ nfc"-i (3.7) 

n=Q ''O n=0 

Our transform I^ is ^„=o n (^n^"'/n\ as before, but if we replace the sharp cut-off 
by a Gaussian one we obtain an alternative to A^ 



Jo U 



Ctfj Tbrh 6 



1 ^ 

-Y^nnmajs- 



(3. 



n=0 



with a transform ^„=o at r(?T-/2) a„A"'/(2n!) which clearly converges much more 
slowly (if at all) due to the T{n/2) in the numerator. That this is the reverse to 



what occurs in the Polchinski flow equation in |3^ is because in that application the 
corresponding expansions involve positive powers of s rather than negative ones so 
that r(ra/2) appears instead in the denominator and enhances convergence. 
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4. Surface Critical Scaling in The Ising Model 

(y9^ theory is the Landau- Ginzburg theory that describes the critical behaviour of 



the Ising model [^, |^. The vacuum wave- functional has a functional integral 
representation as a Euclidean field theory with a boundary corresponding to the 
quantisation surface 



oo 



xl>M= /P0exp( - I dt I rfx(-(^0^ + 0'^ + m^0^j + |i^ 

+ / dx(f o^{x,0)) (4.1) 



where the satisfies Dirichlet boundary conditions 0(x, 0) = 0, and a regulator is 
imposed (which we take as the same momentum cut-off as before), ipo appears as 
a source term coupled to (p so that functional differentiation results in an insertion 
of 0. For a particular value of g the ^ ^ limit of \1/ models the critical behaviour 
of the two-dimensional Ising model with a boundary on which the spins a ~ 
obey Dirichlet boundary conditions, this is the Ordinary Transition. The correla- 
tor (6"(0, 0) (t(x, 0)) can be computed from the exact Onsager solution or from the 
conformal field theory four-point function using the method of images and shown to 
depend on \x\ as l/|a:|, |^. So in the Conformal Field Theory limit of (p^ that part 
of W quadratic in the field,/ dx dinpo{x) (po{y) ^o{x — y) = J dxdy(p{x) ip{y) T{x — y), 
should have F oc l/\x — y\. 

We have studied massive ip'^ theory, but the Renormalisation Group provides a 
means of extracting the massless behaviour |^ , |3^ . It starts from the observation 
that at short distances the functional integral ( [4.1|) expressed as it is in terms of 
bare quantities is approximately independent of rriQ. This is based on a perturbative 
analysis, because in terms of Feynman diagrams differentiating ^ with respect to rriQ 
introduces a 0^ insertion in the bulk involving two propagators that are suppressed at 
large momentum. Now, from the functional integral representation Fq is a function 
of the bare variables go, mo and the cut-off and is related by wave- function renor- 
malisation to F which is given by our solution in terms of fi and g. By dimensional 
analysis Z depends on fi and s only in the combination fis whilst F depends on x — y 
as F{fj,{x — y),g)/\x — yp with F a dimensionless function, so 

Fo(x - y, go, mo, s) = Z'\ij,s, g) F{ii{x - y),g)/\x - y\^ (4.2) 

Since Fq can be taken to be independent of mo for small |x — y| the left-hand-side 
does not change if mo is varied whilst go and s are held fixed. As mg varies /x and g 
will vary together, so that g can be thought of as a function of yU with a /5-function, 
(3 = dg/dlogfi = P{g). Since Fq does not change if we vary /x and g from /xi and gi 
to fj,2 and g2 = 5'(/ii,/i2,5'i)we obtain 

F{i^2{x - y),g2) = 777 ^ F{fii{x - y),gi) (4.3) 
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Now 7 = d log Z/d log fi is a function oi g, so 



— exp / — 7 = exp / —7 (4.4) 



If g,, is a zero of the /5- function and if we assume that 7 is continuous then for gi, g2 



close to 5^* (|4.4|) behaves hke 



exp / 7 X ^// x 7(^J ~ — (4-5) 



92 



(g - 9*) P'ig*) V/^i 



If in addition we assume that F{x, g) is continuous near g^, then the limit as gi -^ g^ 
of f3) is 

F(/i2(x-2/),^,)= r^") F{fi,{x-y),g.) (4.6) 

which implies the homogeneous dependence on x — y: F ~ |x — ^['''(s'*) that is char- 
acteristic of a Conformal Field Theory. 

The correspondence betwen the Ising model and ip^ theory requires that we 
should be able to obtain 7(5'*) = 1 from our solution. The series given in the 
Appendix express go and Z as functions of fi, s and g. Differentiating these with 
respect to log/i whilst keeping go and s fixed gives 

n dgo ^ n dgo dZ dZ 

a log/i og a log/i og 

From which f3 and 7 are obtained as power series in 1/s with the results also given 
in the Appendix up to order 1/s^^. We have calculated these functions for a variety 
of values of g using the transformation of Section 3 to extract the s — > limit. The 
results are shown in Figure ^ and demonstrate that the /9-function has a zero at g 
roughly equal to 2 where 7 is approximately unity in agreement with the Ising model 
result. Studying the region in the vicinity of 5^ = 2 leads to values of 1.955 for the 
zero and 0.961 for the corresponding value of 7. 

Comparison with the exact result shows a departure from it by four per cent, 
but we can estimate the error a priori by examining the shape of I^ . For small A the 
error is mainly composed of exponential terms that are suppressed as A increases, 
we can estimate the dominant term by fitting dl^ /dX to an expression of the form 
exp(mA + c) for < A < 0.9A* (the region is only taken up to 0.9A* rather than A 
because the dl]j /d\ and dl]^ /d\ diverge for smaller A). In Figure |^ we plot I]J and 
exp(mA + c)/m which shows that the exponential terms account for about six per 
cent of the value of /^, this is of course only a rough estimate. 

This approximation of the anomalous dimension depends on the number of terms 
we keep in the series expansions, which in turn determines the type of terms we need 
to retain in our local expansions of W and U . So far we have worked to order 1/s" 
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Figure 7: ijj (upper curve) and exp(?TT,A + c)/m approximating the error. 



with n = 17 which has required including terms with F fields and D derivatives 
such that F + 2D < 40. Since the true continuum limit should be obtained as 
n — >^ oo it is worthwhile to see how the results vary with n. In Figure |^ we plot the 

1-- 
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Figure 8: The approximation to r/ for n = 17, 15, 13, 11, 9 
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approximations obtained by keeping n terms against 1/n for n = 17, 15, 13, 11, 9. 
We also show the exact result as though it would occur at n = oo, to guide the eye. 
The plot shows roughly linear behaviour, and when a straight-line fit through these 
results is extrapolated we obtain a value of 1.003 for the anomalous dimension, which 
is in very good agreement with the exact Ising Model result. 

We can compare this result with the predictions of one-loop perturbation theory. 
The vacuum functional was computed in the semi-classical approximation in |^ 
(the perturbative treatment of the three-dimensional case is given in ^^) where the 



coefficients of ^p^ and ^p'^ in W were found to be -go/{9QM) + 3.973. lO-^^j^/M^ and 
— 1/(4M) + (7o/(487rM), respectively, where M is related to the particle mass by 
M = ft — go/i^S-Kfi). These imply that at one-loop order /3 = —2g + 0.617(7^ and 
7 = g/ii, so that g^, ~ 3.242 and 7(5'*) ~ 1.032, which is slightly more accurate than 
our estimate based on series with 17 terms, but considerably less accurate than the 
extrapolated value of 1.003. 

The expansions in the appendix can be re-ordered to yield estimates for the 
coefficients of the loop expansion. For example, if we retain only those terms up to 
and including second order in g in the beta-function we obtain 



2/I5 1 65 1 , 3773 1 179 



P^ -2g + g'{- ^-r-^ + 



V 8 TT/is 72 fi^s^TT 5760 jJ's^ii 336 ji^s'^n 
341759 1 1112423 1 38151475 1 



746496 /i9s% 2737152 /iiisi% 103514112 ^i^^s^^ti 
182364253 1 9235090201 1 \ 



537477120 ^i^^s^H 29238755328 ^Ji^^ s^H 

Using the Borel resummation technique to find the s — > limit of the coefficient of 
g'^ gives 0.633 which is in reasonable agreement with the one-loop value 0.617 given 
that this series has so few terms in comparison with the previous series we have 
considered. Similarly the one-loop contribution to 7 is approximated by 

/9 1 ,,0 1 63 1 27 1 



8 TT/iS Ijfis^Ti 128 fi^S^TC 64 /i's'TT 

385 1 351 1 10395 1 

+ 



1024 /i9s% 1024 /iiisi% 32768 fx^^s^H 

2431 1 73359 1 \ 

~8192 /iiSgiSTT ^ 262144 fi^'^s^^Tr) ^'^ 

Resummation yields the value 0.336(7 which should be compared with g/n ^ 0.318 (7. 
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We can go further and estimate the two-loop contribution to the beta-function 
from the series 

,/ , 1 65 1 5975 1 1082059 1 

/ -9/2— ^ + T^ 



7r/i2s2 ISvr/iV 10368 7r/i6s6 60480 tt ji^s^ 
637469833 1 73790276321 1 52328568942469 1 

7464960 TT/iiOsio ^ 215550720 ir fi^^s^^ ^ 39852933120 ir fi^^s^^ 
1607134369983293 1 



322808758272 7i fi^^s 



160I6 



(4.10) 



which gives —0.8135'^ on resummation, so to two-loop order the /^-function ceases to 
have a zero. The disappearance (and reappearance) of the zero in successive orders 
of perturbation theory is encountered elsewhere and requires a skilful handling of 
the perturbation expansion, for example with the use of Pade-Borel resummation in 
g, ||3^. This is to be contrasted with the results from the Schrodinger equation for 
which there is always a zero in the /3-function although its position is corrected at 
each increasing order of the approximation. 

5. Conclusions 

We have solved the functional Schrodinger equation for the vacuum of if'^ theory and 
for a state describing a stationary particle by expanding these states in terms of 
local functionals. The solutions were obtained as power series in 1/s, where s is a 
short- distance cut-off. The s ^ limits of these series were found using a variant of 
Borel re-summation that avoids the use of analytic continuation. A standard Renor- 
malisation Group argument was used to compute an anomalous scaling dimension 
in the massless theory. Extrapolating in the number of terms of the series gave the 
result 1.003 for this scaling dimension which agrees well with the exact result of unity 
obtained from the conformal field theory treatment of the surface critical scaling of 
the Ising model at the Ordinary Transition. 

We have only developed the solution to the point where we could make the 
comparison with the Ising model. Further work remains to be done to describe 
the particle dynamics in this approach, for example by checking Lorentz invariance 
and constructing an S-matrix. Neither have we computed inner products of states, 
although their form is determined by the Hermiticity of the Hamiltonian and they 
should also be calculable as re-summed expansions in 1/s. 

We concentrated on two dimensional scalar theory as this allowed the series ex- 
pansions to be obtained using only modest computing resources. Similar resources 
should be sufficient to study other models such as Landau- Ginzburg theories with 
higher order polynomial interactions and integrable models with exponential inter- 
actions where the wealth of exact results would provide useful tests of the method. 
There is nothing intrinsic to the approach that limits it to two dimensions, so it 
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should be possible to apply it to critical phenomena in 3 dimensions and gauge 
theories in three and four. Indeed the method was motivated by the desire to sys- 
tematise a heuristic model of the Yang-Mills vacuum that gave an accurate formula 
for the string tension in three dimensions. Applying it to that case to compute series 
with large numbers of terms would require some effort, so we cannot at the moment 
explain why the formula for the string-tension agrees accurately with Monte-Carlo 
results, despite being based on just the leading term of the local expansion. However 
we can understand why the corresponding calculation in four-dimensions appears to 
give a string-tension that diverges with the cut-off: this is just the first term in a 
power series that should re-sum to a finite result. As in the case of scalar field theory 
standard perturbation theory is obtained by organising the local expansion in powers 
of the coupling. In this context it has already been shown how the usual one-loop 
beta- function for Yang-Mills theory is obtained in this approach ||38||. 
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A. Solutions 

The solutions listed below for the series expansions of Z, qq, j3 and 7 were computed 
up to order 1/ s^'^ using MAPLE, working with rational numbers. The exact values 
for the cofiicients have been evaluated to ten significant figures to save space. 

Solving the Schrodinger equation yields the series expansion of the wave-function 
renormalisation, Z = 'Ylm^n/i^^s)^-, with the first 18 coefficients given by 

zo = 1.0 

zi = -0.1193662073^ 

22 = 0.05066059180/ 

Z3 = 0.03978873576^ - 0.07596748147 / 

Z4 = 0.1658372186/ - 0.04239306467 ^^ 

Z5 = 0.1322116938/ - 0.4702098769/ - 0.02238116387^ 

Z6 = 1.600512277^"^ - 0.4288302602/ + 0.02320184743/ 

Z7 = 0.01492077591 g + 1.643970766 g^ - 6.280656209 g^ - 0.2150394698 / 

zs = -7.118564002/ + 27.72021990/ + 0.03759400389 ^^ + 0.9449811249/ 

zq = -0.01087973244^ - 4.714821300^^ + 34.16960776^^ 

-135.3903272/ + 0.3954218093/ 

zio = 723.5185086^^^ + 25.14159983/ - 179.4495800/ 

-0.2452396545^^ - 2.090535257/ 

z^^ = 0.008392936450 g - 4195.152125 g^^ + 14.06398176 / 

-143.8822880/ + 1021.740867/ -0.9435065158/ 

Z12 = 26222.36179^^2 - 6265.013026 c/^° - 92.19320030/ 

+880.6482723^^ + 0.9593860756 ^^ + 3.790793362/ 

z^3 = -0.006731834444^- 175775.0461 ^^''^ + 41156.19353^^^ 

-45.34168199^^ + 626.1457656^^ - 5746.450353/ 

+3.174250642^^ 

zu = -288442.7166^^2 + 1258090.875^^^ + 39859.60610^^° 

+375.4551915/ - 4438.133512/ - 3.420232517/ 

+9.605144838/ 

2^5 = 0.005556434780 g + 2149122.481 g^^ - 293133. 1749 ^" 
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+138.6466218/ - 3087.933824^^ + 32967.21250^ 

-14.23407683 c/^ - 9578735.713^^^ 

ziQ = 2280136.160^^^ - 16970850.80^^^ - 256999.1253^^° 

-1543.206840 / + 25586.11980 g^ + 11.94349655 / 

-272.9892577/ + 77323626.71 g^^ 

z^j = -0.004688241844^ + 75.31011486^ + 16703.89651/ 

-141.2669974 /+ 141642086.6 ^^^ - 18718142.21 ^^^ 

+2102889.466^^^ - 217344.7165/ - 659842155.7 c/^^ 
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Solving the Schrodinger equation yields the series expansion of the bare-coupling, 

9o = Ai^ J2n i^n/ifJ'S)"', with the first 18 coefficients given by 

Ko = 1-05' 

Ki = 0.1989436788/ 

K2 = -0.03641230036/ 

Kg = 0.08207763547/ - 0.05747261834 c/^ 

K^ = -0.1813715346/ + 0.02286762825/ 

Kg = 0.5276961593/ -0.1344338864/ + 0.02978628969/ 

Kg = -1.827461455 g'^ + 0.4342310123 g^ + 0.02099375315 / 

K7 = 0.2040031164 / - 0.01884175582 g^ - 1.725204731 / + 7.268618386 g^ 

Kg = -0.8376851525/ - 0.1619837103/ + 7.650541972/ - 32.43432692/ 

Kg = -0.3373967928/ + 0.01324798773/ + 4.478764811/ 

-37.43268679/ + 159.8791617^^° 

Kio = 1.190013912/ + 0.6263468208/ - 24.78826072/ 

+ 199.6885243/ - 861.1578312^" 

Ku = 0.7212680080 / - 0.009951254171 g^ - 10.64530019 / 

+ 146.1049618/ - 1152.010150^^° + 5027.712587^^^ 

fi;i2 = 3.152841368/ - 2.174458226/ + 76.02682717/ 

-915.5371834/ + 7143.555824^" - 31617.81483^^^ 

Ki3 = -2.418827270/ + 0.007821150456/ + 13.86178637/ 

-549.6708328 / + 6091.031330 g^^ - 47386.89418 g^^ 

+213090.6078 c/^^ 

Ku = -68.20413916/ + 7.396453034/ - 178.3373919/ 

+4076.316308/ - 42945.16713^" + 334961.5945^^^ 

-1532576.001 c/^^ 

Ki5 = 11.95109263/ - 0.006353030446/ + 151.2121465/ 

+ 1867.345135/ - 31358.18697^^° + 320273.7334 ^^^ 

-2514689.736^^^+11719597.0^^^ 

^16 = 819.5644108^^ - 25.22638048/ - 790.3286415 c/^ 
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-17670.92565 g^ + 251437.2784 g^^ - 2521699.943 g^^ 

+ 19992199.42^^^ - 94979792.35^^^ 

Kn = -167874034.0 c/^^ - 2682.130056/ + 0.005291499834 c/^ 

-69.91144120 / + 436.3776317/ + 163982.6175 g^^- 

2105884.503^^2 + 20923210.23^^^ + 813423914.7^^^ 
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The solutions for Z and g^ lead to the series expansion of the /3-function, /3 
X]„/5n/(/Us)", with the first 18 coefficients given by 

/9o = -2.0(7 

/9i = 0.5968310364/ 

/52 = -0.4559453262 c/2 

/Sg = 0.9853095730/ - 0.2873630916 ^^ 

/54 = 0.3658820519/ - 2.814252613/ 

/^s = 0.2085040278 g^ + 9.833762561 ^^ - 1.749096318 / 

/Se = 0.05839063195^^ - 39.74109388/ + 7.374962709/ 

/57 = -0. 1695758024 /+ 180.4235628/ -34.76674176/ + 2. 702926953/ 

/5g = -903.9912677/ - 1.812756260^^ + 178.4867048^^ - 15.43817775^^ 

/39 = 0.1457278651/ -990.1865501/ + 96.47016296/ 

-4.096644239 / + 4940.346009 g"^ 

/5io = 5898.620256^^ - 29206.19354^" + 8.652316689/ 

-615.2794259 g'' + 24.26339449 / 

/5ii = -0.1293663043 g^ + 4091.802192 / - 241.2179221 / 

+7.386237270/ - 37554.52475^^° + 185620.8411 ^^^ 
/5i2 = -28511.70908/ + 254595.6560^" - 34.68565605/ 

+ 1992.597606^^ + 59.07548891/- 1262105.326/^ 

/5i3 = 9144064.399/^ + 0.1173172568/ - 16205.08595/ 

+289.2141 186 / - 26.54957599 / + 208570.6632 c/^° 

-1832311.357/2 

/3i4 = 133241.9688/ - 1602568.994^" + 133.0389542 ^^ 

-4843.070749^^ - 1489.424050^^ + 13962796.56^^^ 

-70350835.71^^^ 

/5i5 = -112396967.2^1^ - 0.1080015175/ + 57364.96045/ 

+4603.125796 / + 573028054.8 g^^ + 177.9010353 / 

-1123852.943^1° + 12930212.94 ^^^ 

/5i6 = -602711.1263/ + 9794809.468^" - 504.4372323/ 
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-24799.33347^^+19590.50509/ - 109478969.7^^^ 

+953696332.6^^^ - 4928124568.0^^^ 

(3n = 44638767880.0^^^ + 6126484.834 ^^^ - 82996.74004 c/^ 

-1336.249746/ + 0.1005384968/ + 971831675.9^^^ 

+27515.67004/ - 88542747.03^^^ - 8512571797.0 ^^"^ 



-27- 



The solutions for Z and (70 lead to the series expansion of the the logarithmic 
derivative of Z, 7 = X]n7ra/(/^s)", with the first 17 coefficients given by 

71 = 0.3580986219^ 
72 = -0.3324601337^2 

73 = -0.1989436788^ + 0.7407774312^^ 

74 = -2.157405193^ + 0.3591976683/ 

75 = -1.514499301^ + 7.609034060^ + 0.1566681471^ 

7g = -30.90102012/ + 6.324551298/ - 0.2470671327/ 

77 = 2.810587640 / + 140.6420098 g^ - 29.47246351 / 

-0.1342869832^ 

78 = -15.56377334/ - 0.4393122678/ - 705.4720543/ 

+ 149.9640642/ 

79 = 3856.625518 / + 0.1196770568 g - 5.755944500 / 

-825.7869756 g^ + 92.59293353 g^ 

710 = 3.423764699 g^ + 4887.996878 g^ - 573.0603460 g^ 

+38.04564185 / - 22794.88025 g^^ 

711 = -0. 1091081738 g - 30947.00771 / + 144798.2478 g"^ 

+15.15984462 / + 3727.952024 g^ - 299.6137729 g^ 

712 = 208764.5292^^° - 983835.0583 ^^^ - 15.34219367/ 
-25536.87550^^ + 2257.479359^^ - 75.61099893/ 

713 = 0.1009775167^+184253.1860/ - 1495814.945^^^ 
-56.87562475/ - 17294.46868/ + 1035.456051/ 

+7122090.122^^^ 

714 = 61.55740610/ + 136372.4752/ - 9808.072105/ 

-54746191.69^^^ - 213.8081010/ - 1399600.487^^° 
+ 11353026.88^^2 

715 = -0.09445939125 g - 3285.663980 g^ - 91057448.51 g^^ 

+445520492.1^^^-1114679.686/+ 11183177.87^^^ 

+286.2007410 / + 90432.28371 g^ 
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716 = 42485.83958 / + 6557.967853 g* + 9477782.019 g^^ 

-93893690.17^^2 - 3828085628.0^^^ + 770074284.4^^^ 

-238.8640816/ - 829515.1962/ 

717 = -514466.9761 g'^ - 83973087.42 ^" + 0.08907659505 g 

-1685.504320^^ + 1614.901156^^ + 827365423.6^^^ 

-6852797795.0 g^^ + 7721868.666 / + 34643996560.0 g^'^ 
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